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Abstract 

We test the compatibility of chiral symmetry, dynamic mass generation of the nu- 
cleoli due to spontaneous breaking of chiral symmetry, and the description of finite 
nuclear systems by employing an NJL model understood as a chiral invariant effec- 
tive theory for nucleons. We apply the model to nuclear matter as well as to finite 
nuclei. In the latter case, the model is adjusted to nuclear ground-state observables. 
We treat the case of a pure chiral theory and the physically more realistic case where 
a portion of the nucleon mass (160 MeV) explicitly breaks chiral symmetry. The best 
version of this current model is found to deliver reasonably good results simultane- 
ously for both finite nuclei and the nucleon mass, which supports our motivation of 
probing a link between low-momentum QCD and the nuclear many-body problem. 
However, the observables calculated for finite nuclei are not as good as those coming 
from existing relativistic mean field models without explicit chiral symmetry. 

Key words: Nambu-Jona-Lasinio model, relativistic mean-field model, nuclear 

matter, finite nuclei 
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1 Introduction 



QCD is believed to be the fundamental theory of the strong interaction. It 
has non-perturbative behaviour at low energies. One of its important non- 
perturbative phenomena is the nontrivial vacuum structure: the vacuum is 
populated by scalar quark-antiquark pairs leading to a finite expectation value 
of (qq), the chiral or quark condensate. The Nambu-Jona-Lasinio (NJL) model 
[1,2,4,3], originally invented for nucleonic, but mostly used for quark degrees 
of freedom, schematically implements this rearrangement of the vacuum. The 
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quark mass is generated by spontaneous symmetry breaking due to the neg- 
ative energy states of the Dirac sea which contribute to the scalar density. 
In the context of nucleonic degrees of freedom, the major part of the nucleon 
mass is generated by the chiral condensate, which is modeled (see below for 
a more specific explanation) by (NN), where N denotes nucleon states. A 
smaller part of the nucleon mass emerges from an explicit chiral symmetry 
breaking contribution. 

A link between hadronic effective field theories for nuclear matter and finite 
nuclei and QCD can be established via the symmetries of (massless) QCD: 
predominantly Lorentz invariance and chiral symmetry. There exist models 
employing a nonlinear realization of chiral symmetry that have proven to be 
successful [25,12]. Successfull models for nuclei not considering explicit chi- 
ral invariance are relativistic mean-field (RMF) models, either the Walecka- 
type (RMF-FR) with explicit finite range meson exchange [5,6] or the point- 
coupling type (RMF-PC) with contact interactions and derivative terms to 
simulate finite range [7,8,10]. Both of these deliver accurate descriptions of fi- 
nite nuclei observables such as binding energies, densities, radii, surface thick- 
nesses, and ground-state deformations. They can be considered compatible 
with chiral symmetry in the sense that power counting (NDA - naive dimen- 
sional analysis) based upon chiral symmetry works. This is more transparent 
in the point-coupling variant of these models. Thus, nuclei are compatible with 
QCD scales [13,8,14]. 

By describing nuclear matter with NJL models one employs effectively (and 
schematically) three features of QCD, namely, chiral symmetry, the mass gen- 
eration due to its explicit and spontaneous breaking, as well as the binding 
of nuclear matter (i.e. the nucleon in-medium self-energy) resulting from the 
change of the condensate. Nonlinear NJL models for nucleon degrees of free- 
dom have been applied to nuclear matter [16,22,23,24] and, in a simple version, 
to finite nuclei [17]. Nucleonic NJL models provide an attempt to incorporate 
chiral symmetry and its dynamical breaking within the framework of relativis- 
tic mean-field models. The chiral condensate is modeled in terms of nucleonic 
degrees of freedom, (qq) — ■> (NN). Some comments on this relation are in or- 
der. It has been shown by several authors [25,26,27,28,29] that expressing the 
quark Dirac sea by the nucleon Dirac sea is inapproriate. Strictly speaking, 
the nucleonic degrees of freedom, which are well suited for the valence par- 
ticles in the spirit of effective field theory, are the wrong ones for describing 
the chiral condensate. Replacing (qq) by (NN) is not compatible with the 
large N c limit of QCD. Furthermore, the composite nature of the nucleon sup- 
presses corresponding vacuum loops. In our case, we do not think of a literal 
replacement. The sole model property needed stemming from (NN) is a large 
(negative) contribution to the scalar density (the scalar mean field) which in 
our best-fit model contributes to about 80 % of the nucleon mass. This usage 
of the nucleon Dirac sea in finite systems is put to the test. 
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Based on our experience with RMF-PC models [13,8], we would like to test an 
extended and physically more realistic version for finite nuclei, including im- 
portant ingredients such as derivative terms, the Coulomb force, and the center 
of mass correction. Nuclei provide a more stringent test of this model than in- 
finite nuclear matter does since the shell structure and density distribution are 
quite sensitive to the size and the density dependence of the mean-field poten- 
tials. Our basis for comparison will be the RMF-PC force PC-F1 determined 
in Ref. [8] which respresents the accuracy that can be obtained in modern 
mean-field approaches. We call the model used here, to distinguish it from the 
NJL model for quark degrees of freedom, the nucleon NJL (NNJL) model (we 
correspondingly denote that for quarks by qNJL). Note that a model for the 
deuteron has also been denoted as NNJL [9]. 



An NJL approach for nuclei resembles the RMF-PC model in the Hartree ap- 
proximation which we also use here, since all densities and currents containing 
75 matrices vanish in ground states of even-even nuclei. The most prominent 
difference between the two is that RMF models contain an explicit mass term 
and are treated in the no-sea approximation, whereas the mass is dynamically 
generated within NJL models. Another difference is that in contrast to the 
linear RMF-PC model, the linear NJL model cannot describe saturation of 
nuclear matter (if reproducing the vacuum nucleon mass at the same time). 
Saturation can only be obtained at the price of a small or vanishing nucleon 
mass [15]. Nonlinear NJL models, however, are capable of simultaneously de- 
scribing saturation and the vacuum nucleon mass, as demonstrated by Koch et 
al. [16], as well as the chiral phase transition. Another important difference is 
that imposing chiral symmetry in the NJL model lagrangian restricts allowed 
terms. The building blocks of the NJL model that in the Hartree mean-field 
approximation reduce to isoscalar terms are the chiral invariant combinations 



(ipip) 2 - (ip-y 5 Tip) 2 and {ipi^f + (^75^7/^) 



This constrains the density 



dependence of the effective mass (and hence the scalar potential). 



With our approach, which to our knowledge is the first realistic NJL model 
for light and heavy finite nuclear systems in the mean-field approximation, 
we wish to probe a linkage between low-momentum QCD and the properties 
of finite nuclei by simultaneously employing a chirally symmetric interaction 
lagrangian and mass generation by spontaneous breaking of chiral symmetry. 
We do not intend to create an effective field theory for nuclei that is better 
than the already established models, like the relativistic mean-field model and 
its chiral versions. 
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2 The model 



The starting point of our investigations is based upon the nonlinear chiral 
invariant model of Koch et al. [16,17], which we formulate in its SU(2)y <8> 
SU(2)a form (where mo = 0): 



t-'K — £free + £ 



-4/ 



c 



8/ 



•-/ree 



A/ = gi 

A/ = 93 



{4)ip) 2 - (^ripf 



~92 



X 



(2) 



The cross-density term associated with the coupling constant g 3 is usually 
missing in RMF models since saturation is achieved already in the linear 
version. Not so in the NNJL model: here the cross term is necessary to achieve 
saturation as well as a realistic vacuum mass at the same time (for a discussion 
of the problem of matter stability within the NJL model see [30]). In the 
original version from Koch et al., where mo is set to zero, the model has no 
explicit symmetry breaking terms. In our analysis, we consider the realistic 
case with an explicit symmetry breaking term (m > 0). For comparisons 
to the model of Koch et al. in nuclear matter we also consider the limit of 
massless nucleons. 



We extend Lk for nuclear matter with the following terms: 



£exti — 9 a 



95 



(3) 



This is motivated by the fact that the effective mass in the model of Koch et 
al. [16] is too large (m*/m vac = 0.92) and must be smaller to generate spin- 
orbit splittings in nuclei of the right magnitude. In contrast to the RMF-PC 
approach, the parameters gi, g^, and g± not only govern the scalar potential 
felt by the nucleons, but simultaneously contribute to the nucleon mass at 
zero baryon density. This is an important additional constraint which does 
not appear in the case of RMF-PC. 

For finite nuclei, in a similar fashion, an isovector-vector term is added (hav- 
ing the quantum numbers like the p-meson) as well as derivative terms. The 
derivatives impose no problem since chiral symmetry is a global symmetry. 
They are an important and necessary ingredient in RMF-PC approaches and 
similarly in the Skyrme-Hartree-Fock (SHF) approach [31] to simulate finite 
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range and to characterize the surface region of the nucleus. We also add the 
Coulomb field. The extension reads: 



(4) 



When applying to finite nuclei, these terms are necessary for RMF-PC models 
and constitute a minimal extension as shown in [7,8]. The parameters of the 
model are the nine coupling constants gi, g 2 , 93, #4, g$, ge, #7, gs, gg as well as 
the (sharp) cutoff A c (see below). We regard this approach as a phenomeno- 
logical model in which all coupling constants are determined by least-squares 
adjustments to measured nuclear ground-state observables. We do not take 
into acount low-energy theorems. One could, in principle, derive a generalized 
Gell-Mann-Oakes-Renner (GOR) relation as can be done in linear g—uj models 
[18] by assuming that the vacuum expectation value of the explicitly symmetry 
breaking term is equal to the corresponding term in low- momentum QCD, see 
also [15,24]. The Goldstone boson in our model, however, is an antinucleon- 
nucleon state and thus can hardly be identified with the pion. We feel that 
these questions (at least at this stage) ask too much of this phenomenologically 
based approach. 

The NNJL model resembles the phenomenological linear a—uo model [18,19,20,21] 
when replacing the fields with their lowest order approximation, namely a — > 
ipif} and 7r — > 'ijj^Txjj (in finite systems additional derivative terms become nec- 
essary). The spontaneous breaking of chiral symmetry in linear a — uo models 
is usually realized by assigning a finite expectation value to the a field due to 
a mexican-hat type potential and then shifting the field. In our case, the NJL 
mechanism leads to a finite vacuum expectation value of the scalar density. 
Both models share the linear realization of chiral symmetry. 



3 Nuclear matter 

Our nuclear matter model is 



nm 



(5) 
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containing 5 coupling constants. Derivative terms as well as isovector terms 
do not contribute in infinite symmetric nuclear matter, and nucleons are mo- 
mentum eigenstates (plane waves). The effective mass is given by 

m* = m - 2(#i + g 3 p 2 v + 2g 4 p 2 s )p s (6) 

In this equation, p v denotes the vector density and ps denotes the scalar 
density, consisting of the valence part and the negative energy states: 



Ps 



f kF *u m * _ _J_ f Ac d * k m * ( 7 ) 
(27r) 3 io y/k 2 + m* 2 (27r) 3 io y/k 2 + m* 2 



where 7 is the spin-isospin degeneracy factor (7 = 4 in symmetric nuclear 
matter) and kp is the Fermi momentum of the valence (positive energy) nu- 
cleons. We can imagine the sharp three-momentum cutoff A c as separating 
the active negative energy states from the ones that do not contribute to the 
scalar density. Stated differently, the effective interaction is only sensitive to 
states up to the momentum A c . We call this the no deep-sea approximation. 

The binding energy per particle B/A is given by 

B/A = (Too - TSS^j /p v - m vac (8) 

where T o is the 00 component of the canonical energy-momentum tensor. Tqq c 
denotes the vacuum contribution that can be obtained by setting the vector 
density to zero and using the vacuum nucleon mass m vac . We obtain 

rk F „ / rA c 



Too = / F d 3 kJk 2 + m* 2 - f " d 3 k\/k 2 + m* 2 
Jo Jo 

+ gips 2 + g2pv 2 + gzps 2 pv 2 

+ + g5pv A (9) 



and 



^ = " lo C + ^ + + ^Svac ( 10 ) 

Note that the constant shift of the negative energy states caused by the vector 
potential does not contribute since it is becomes zero by the normal ordering 
of the baryon number operator. The integrals for the scalar density and the 
two binding energy contributions are done analytically. The energy computed 
within the NNJL model results from two large cancellations: one is, as in the 
RMF approach, the cancellation of large attractive scalar and repulsive vector 
fields and their contribution to the total energy, while the other is the difference 
of the contribution at finite baryon density and the vacuum contribution from 
the active Dirac sea states. 
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From the experience with RMF models we know that one of the attractive 
features of this approach is that the spin-orbit splittings of nuclei are nicely 
reproduced, courtesy of the large (several hundred MeV) scalar and vector 
potentials. Since with C K) the effective mass is about m* /m vac = 0.92 at sat- 
uration [16], the spin-orbit splittings in finite nuclei can be expected to be 
far too small (we will demonstrate this statement in the next section). To 
get closer to a model which also has a satisfactory predictive power for fi- 
nite nuclei, especially their shell structure, several least-squares adjustments 
were performed to pseudo-observables p v = 0.17 fm -3 (0.5), E/A = -16.5 
MeV (1.0), m vac = 938.9 MeV (1.0), and m*/m = 0.6 (2.0) [the numbers in 
parentheses indicate the uncertainties in percent entering the x 2 ]- The uncer- 
tainties are, of course, rather arbitrary and reflect only the accuracy that is 
expected of the pseudo-observables. The baryon density p v = 0.17 fm~ 3 has 
been favoured over p v = 0.16 fm -3 since this value has been used in the adjust- 
ments of Ref. [16]. For the adjustment N-l a compressibility of K = 240 MeV, 
with 5 % error, has been included in x 2 ■ 

Our results are summarized in Tables 1 and 2. These adjustments have been 
performed with various fixed values for A c since it turned out that when in- 
cluding it as a free parameter it does not experience significant changes. Fur- 
thermore, an explicit chiral symmetry breaking mass term improved the x 2 i n 
the adjustments of N-l and N-2 compared to N-K*. According to [32], where 
the deep-inelastic momentum sum rule and the trace anomaly of the QCD 
energy-momentum tensor have been used to separate different portions of the 
nucleon mass, this mass term should not be larger than m = 160 MeV. It 
turned out, however, that larger values are more favorable than smaller ones, 
so the adjustments N-l and N-2 have mo set to 160 MeV. This is a first trial 
value that has also been employed in Ref. [15]. A different choice would be the 
pion-nucleon sigma term ( X ps 50 MeV) which we postpone to forthcoming 
work. 

Table 1 

Bulk properties of nuclear matter for the set N-K, the parameter sets obtained in 
least-squares adjustments to nuclear matter properties (N-K*, N-l, N-2), and to 
finite nuclei (N-FN), as well as for the RMF-PC force PC-F1. 



Set 


m vac [MeV] 


PV [fm" 3 ] 


E/A [MeV] 


m*/m vac 


K [MeV] 


N-K 


942.4 


0.172 


-16.6 


0.92 


354 


N-K* 


935.7 


0.167 


-16.6 


0.82 


355 


N-l 


937.1 


0.172 


-16.3 


0.66 


311 


N-2 


935.0 


0.178 


-16.5 


0.72 


349 


N-FN 


941.2 


0.155 


-16.4 


0.75 


388 


PC-F1 


938.9 


0.151 


-16.2 


0.61 


270 
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Table 2 

Coupling constants of the parameter sets discussed in the text. 



oet 


rp 2l 


92 [fm z J 


rp 8i 


rp 8i 
94 [im 5 J 


rp 8i 
95 inn 5 ] 


a rc — li 
A c [tm L \ 


mo [MeVJ 


N-K 


1.16 


3.40 


0.92 


0.0 


0.0 


3.25 


0.0 


N-K* 


1.5916 


3.1950 


1.3665 


0.0773 


0.0 


2.75 


0.0 


N-l 


3.3064 


3.9531 


2.2270 


0.9910 


-0.8917 


2.00 


160.0 


N-2 


1.7678 


3.1279 


1.0908 


0.4685 


0.0 


2.35 


160.0 




p v [fm" 3 ] E [MeV/fm 3 ] 

Fig. 1. Energy per particle (bottom left) and effective mass (top left) as a function 
of nucleon density as well as pressure vs. density (right) for the six models discussed 
in the text. The causal limit is given by the line P=E. 

Clearly, an explicit mass term and a cutoff below 2.40 fm -1 (N-l and N-2) 
achieve an effective mass which is close (but still larger) to those determined 
in RMF models (m* is larger with tuq = 0, see N-K*). But by comparing N-K* 
to N-K we notice that the additional terms and the inclusion of the effective 
mass in the adjustment procedure drive the model into the right direction 
already. The new sets N-l, N-2, and N-K* tend to a smaller vacuum mass, a 
higher saturation density and a rather high compressibility compared to the 
desired values. The energy at saturation appears to be reproduced by any of 
the sets of coupling constants. It is curious that the model delivers a high 
compressibility regardless of the number and/or values of the parameters. A 
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lower compressibilitiy, K=311 MeV, could only be achieved with a very low 
cutoff value of A c = 2.0 fm" 1 in the set N-l. Overall, N-l performs best due 
to the explicit mass term, the coupling constant g 5 and the inclusion of K in 
the adjustment procedure. The NNJL model appears to favor cutoff values 
of approx. 2.0-2.35 fin -1 . This is still larger than the fermi momentum of 
the valence nucleons, but clearly smaller than values used in qNJL models 
3 fin -1 ). This empirical fact, which has been found also by other authors 
[23,24], still lacks a clear explanation. The cutoff, however, is of the order of 
the a meson mass used in RMF models 450 — 550 MeV). 

Figure 1 shows the effective mass and the energy/particle as functions of the 
baryon density as well as pressure versus energy. Clearly, models with (PC- 
Fl, N-l, N-2) and without (N-K, N-K*) an explicit mass term exhibit a quite 
different dependence of the nucleon mass with density. In a pure chirally sym- 
metric theory, there is a second order phase transition from massive nucleons 
to massless particles. However, in the physically realistic models with m > 0, 
a smooth transition takes place which resembles the behaviour of the effective 
mass in RMF models. For a successful model of nuclei in terms of spectral prop- 
erties, the effective mass needs i) m*/m ps 0.6 and ii) a density dependence 
that leads to the proper radial dependence of the spin-orbit potential [which 
is proportional to the spatial derivative of Vg — VV, where Vg (Vy) is the scalar 
(vector) potential]. We note that the RMF-PC force (PC-F1) and the NNJL 
sets adjusted with respect to nuclear matter pseudo observables (N-K, N-K*, 
N-l, N-2) obey causality. Causality is violated above E >^ 1500 MeV/fm 3 , 
however, for the set N-FN which has been adjusted to nuclear ground-state 
observables (see the discussion in the next section). 

We also note from the lower portion of Fig. 1 that a negative value of g$ 
(set N-l) softens the equation of state due to a reduction of the repulsion by 
the vector potential. At high densities, however, it leads to an instability of 
the equation of state bending downwards again. This is a common feature of 
PC-F1 and N-l, and occurs at such high densities (pv > 1.0 fm -3 ) that it 
does not limit the application of these models to ground-state observables of 
finite nuclei with densities (pv ~ 0.17 fm~ 3 ) that are very near and below the 
saturation density in the equation of state. 



4 Finite nuclei 



Our model for finite nuclei is 

£fn — £>K + C-extl + £ext2 (H) 
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containing 9 coupling constants. Finite nuclei are calculated in coordinate 
space employing matrix multiplications as derivatives in Fourier space. The 
solution with lowest energy is determined using the damped gradient step 
method [33]. Pairing for the nucleons is not yet invoked. The negative en- 
ergy states are, as in [17] and for the previous nuclear matter case, treated in 
the LDA approximation with a spatially dependent effective mass m*(r). The 
positive energy states are calculated explicitly and selfconsistently. The pa- 
rameters containing isovector terms and derivative terms (both isoscalar and 
isovector) do not contribute to the Dirac sea states since they vanish in sym- 
metric homogenous nuclear matter. We employ the microscopic center of mass 
(cm.) correction i? c . m . = < ^ n ^ ^ or ^ ne positive energy nucleons only. Since the 
negative energy states are locally treated as chunks of nuclear matter, their 
translational invariance renders a zero cm. correction. 

Starting from the parameter set N-2 we performed least-squares adjustments 
to ground-state observables of 7 nuclei and to the vacuum mass m vac of the 
nucleon simultaneously. The resultant parameter set we label N-FN. The four 
observables used were binding energy E (uncertainty 0.2 percent), diffraction 
radius R (uncertainty 0.5 percent), rms radius r (uncertainty 0.5 percent), 
and surface thickness a (uncertainty 1.5 percent). The nuclei and observables 
chosen were 16 (E, R, a), 40 Ca (all four), 48 Ca (all four), 88 Sr (E, R, r), 
90 Zr (all four), 132 Sn (E) and 208 Pb (all four). The uncertainty for the vacuum 
nucleon mass was chosen to be 1 % (^ 10 MeV). Using a much larger or 
smaller uncertainty affected the adjustment only to a minor extent - the mass, 
in conjunction with nuclear ground-state observables, appears to be a robust 
prediction. This set of observables (except for the nucleon mass) is a subset 
of the ones used to extract the RMF-PC force PC-F1. 
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Table 3 

Observables of the least-squares adjustment procedure: binding energy (E), diffrac- 
tion radius (R), rms radius (r) and surface thickness (a) (units: [fm] for the radii and 
surface thickness, [MeV] for the energy). The last line indicates the total \ 2 with 
respect to this set of observables. Note that N-K and PC-F1 have been adjusted 
earlier to other sets of observables. 



observable 


expt. 


N-K 


N-FN 


PC-F1 




938.9 


942.4 


941.2 


938.9 


16 0: E 


-127.6 


-191.1 


-127.3 


-127.7 


16 0: R 


2.777 


2.646 


2.805 


2.769 


16 0: a 


0.839 


0.601 


0.790 


0.850 


40 Ca: E 


-342.1 


-445.0 


-345.0 


-344.9 


40 Ca: R 


3.845 


3.693 


3.862 


3.839 


40 Ca: r 


3.478 


3.173 


3.418 


3.453 


40 Ca: a 


0.978 


0.706 


0.893 


0.967 


48 Ca: E 


-416.0 


-530.5 


-412.1 


-416.0 


48 Ca: R 


3.964 


3.833 


3.945 


3.945 


48 Ca: r 


3.479 


3.252 


3.434 


3.444 


48 Ca: a 


0.881 


0.712 


0.867 


0.885 


88 Sr: E 


-768.5 


-932.1 


-765.2 


-769.0 


88 Sr: R 


4.994 


4.861 


4.982 


5.005 


88 Sr . r 


4.224 


3.985 


4.168 


4.197 


90 Zr: E 


-783.9 


-941.6 


-779.9 


-785.3 


90 Zr: R 


5.040 


4.835 


4.997 


5.025 


90 Zr . r 


4.270 


4.012 


4.22 


4.245 


90 Zr: a 


0.957 


0.796 


0.939 


0.943 


132 Sn: E 


-1102.9 


-1376.5 


-1107.1 


-1102.8 


208p b . E 


-1636.4 


-2004.4 


-1642.5 


-1636.9 


208 Pb: R 


6.806 


6.775 


6.742 


6.808 


208p b . r 


5.504 


5.330 


5.438 


5.501 


208p b . a 


0.900 


0.341 


0.852 


0.876 


x 2 




156728 


174 


35 
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3 4 
r[fm] 



Fig. 2. Charge densities of the nuclei 16 (left), 48 Ca (middle) and 208 Pb (right) 
with the forces as indicated. The experimental data are from Ref. [36]. 




Fig. 3. Proton (left) and neutron (right) single-particle levels in 16 0. On the left the 
experimental values for the p-states are shown. Forces from left to right are PC-F1, 
N-FN, and N-K. 



Table 3 lists the experimental and predicted values of these observables for the 
sets of coupling constants fitted to finite nuclei as well as for the set determined 
by Koch et al [16]. Also, the total x 2 f° r the different sets is displayed. We 
can see that the parametrizations N-FN constitute a major improvement over 
N-K due to the additional terms and the inclusion of ground-state observables 
in the adjustment of the coupling constants. The original N-K model delivers 
radii and surface thicknesses that are too small. Due to the derivative and 
additional higher order terms, these observables are described better with the 



12 



Table 4 



Sets of (rounded) coupling constants, cutoffs and explicit masses used in this study. 
The fourth and fifth columns contain the naturalized coupling constants for two 
values of the QCD large-mass scale A. 





N-K 


N-FN 


N-FN (A = 770 MeV) 


N-FN (A = 464 MeV) 


9i 


1.16 [fm 2 ] 


1.775 [fm 2 ] 


0.390 


0.390 


92 


3.40 [fm 2 ] 


3.118 [fm 2 ] 


0.685 


0.685 


93 


0.92 [fm 8 ] 


1.170 [fm 8 ] 


0.189 


0.066 


94 


0.0 [fm 8 ] 


0.471 [fm 8 ] 


0.076 


0.028 


95 


0.0 [fm 8 ] 


2.787 [fm 8 ] 


0.450 


0.163 


96 


0.0 [fm 4 ] 


-0.123 [fm 4 ] 


-0.411 


-0.149 


97 


0.0 [fm 4 ] 


-0.138 [fm 4 ] 


-0.462 


-0.168 


98 


0.0 [fm 2 ] 


0.562 [fm 2 ] 


0.494 


0.494 


99 


0.0 [fm 4 ] 


0.167 [fm 4 ] 


2.234 


0.811 


A c [fm- 1 ] 


3.25 


2.35 


2.35 


2.35 


m [MeV] 


0.0 


160.0 


160.0 


160.0 



new NNJL force N-FN. The force N-K delivers a large amount of overbinding 
which is also cured in the new adjustment. Nevertheless, N-K, which has been 
adjusted to nuclear matter properties only, is able to qualitatively describe 
finite nuclei. The average error with N-FN for most observables lies in the 
range of a few percent or below. The surface thicknesses pose the greatest 
problem while the nucleon mass, which is a prediction and thus an additional 
demand on the model, is reproduced to within 0.24 %. 

We see, however, that the RMF-PC force PC-F1 is superior (its x 2 is smaller 
by a factor of 5 than that of N-FN) for the ground-state observables. Pairing 
(which is included in PC-F1) would enhance the results from N-FN, but not 
dramatically. The major deficencies of N-FN are the too large effective mass 
and compressibility. 

The charge densities of the nuclei 16 0, 48 Ca and 208 Pb are shown in Fig. 2. 
The large density fluctuations of N-K have been cured in the new force, N-FN, 
though they are still too large in lead. Also, the surface region that is much 
too small in a model without derivative terms, like N-K, is predicted to be 
larger with N-FN, though relativistic models still generally underestimate this 
observable [8]. PC-F1 comes closest to the experimental data, N-FN performs 
second best. 

Fig. 3 shows the proton and neutron single-particle levels of 16 O for the forces 
discussed as well as the experimental values for the p states. As has been 
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inferred from the nuclear matter properties of N-K, the spin-orbit splittings 
with N-K come out way too small (the p-levels are almost degenerate), while 
they have the right size with PC-F1. N-FN performs second best with the 
splitting being too small by approximately a factor of two. This relates to the 
effective mass (0.75 vs. 0.61 for PC-F1) that is still a bit too large to deliver 
the correct spin-orbit potential. Thus, the NNJL model displays deficiencies 
similar to the linear a models [20,21] The linear realization of chiral symmetry 
and its breaking by the NJL mechanism (as well as the prediction of the 
nucleon mass) constitute strong constraints on the denstity dependence and 
size of the scalar and vector potentials which seem to prohibit a degree of 
accuracy which has been reached with modern relativistic mean-field models. 
Chiral models with a nonlinear realization [25,12], however, deliver an accuracy 
comparable to modern RMF approaches. 

Table 4 shows the coupling constants and other parameters of the two models. 
For the set N-FN we note that the isovector coupling constant gg has approx- 
imately the value of that of PC-F1, which is g$ = 0.675. This is due to the 
fact that g 8 is sensitive only to the difference of neutron and proton vector 
densities and is thus unaffected by the size of the scalar density. The other 
4-fermion coupling constants have different values due to the fact that the 
scalar density is much different from that of the RMF-PC approach given the 
inclusion of the condensate in terms of nucleonic degrees of freedom. 

As stated in the Introduction, QCD scaling provides a valuable test of the 
mean-field coupling constants and hence of the compatibility of a model with 
low-momentum QCD. It has been applied successfully in nuclear physics (see 
[13] for the first application). Naive dimensional analysis is used to obtain 
dimensionless coupling constants. We use the scaling procedure of Manohar 
and Georgi [34] but without pion fields (they vanish in the mean-field approx- 
imation), which reads 



C ~ -C tn 



Iplp 



/. 2 A 2 (12) 



where £ is a generic lagrangian term in the Weinberg expansion [39]. Here 
■ip are nucleon fields f n and m n are the pion decay constant, 92.5 MeV, and 
pion mass, 139.6 MeV, respectively, A = 770 MeV is the QCD large-mass 
scale taken as the p meson mass, and are usual derivatives. Dirac matrices 
and isospin operators (we use t here rather than r) have been ignored. Chiral 
symmetry demands [35] 



A = / + n-2>0 (13) 
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such that the series contains only positive powers of (1/A). A natural la- 
grangian [34] should lead to dimensionless coefficients q„ of order unity. Our 
more stringent definition [8] is that a set of QCD-scaled coupling constants 
is natural if their absolute values are distributed about the value 1 and the 
ratio of the absolute maximum value to the absolute minimum value is less 
than 10. Thus, all information on scales ultimately resides in the Q mn . If they 
are natural, QCD scaling works. In PC-F1, which has been obtained in an 
unconstrained optimization procedure with respect to nuclear ground-state 
observables, the 9 coupling constants are all natural in terms of our stringent 
definition. 



We have scaled the coupling constants of N-FN in two different ways. Firstly, 
we used the QCD mass scale of A = 770 MeV (the p meson mass) which has 
been used to scale the coupling constants of PC-F1. Secondly, we scaled them 
with the scale A = 464 MeV= 2.35 fm _1 which corresponds to the momen- 
tum cutoff of our model and thus constitutes a naturally occuring scale. The 
resulting scaled coupling constants are shown in Table 4. Note that the scaled 
4-fermion coupling constants are identical in both procedures since the mass 
scale does not explicitly enter in their scaling (this is not true for the higher 
order terms). In both cases, the total set of coupling constants does not fullfill 
our criterion of naturalness. This is due to the ratio of the scaled values of 
#4 and g 9 which exceeds 10. Additionally, the absolute values do not scat- 
ter around 1 but are overall smaller. The scaled constants are approximately 
uniform and are, however, almost natural, that is, the deviations from 1 are 
approximately uniform and are not large. We attribute their suppression to 
an incorrect magnitude and density dependence of the potentials which also 
yields an effective mass that is too large. This result shows that in the NNJL 
model QCD scales are not realized as well as in the established RMF mod- 
els, which is probably connected to the way chiral symmetry is both realized 
and broken here (RMF models can be viewed as consistent with a nonlinear 
realization). 



We note that with N-FN a fair description of nuclear ground-state observables 
and simultanesously the nucleon mass occurs even with a high value of the 
compressibilty (K=388 MeV). The nuclear ground-state observables do not 
appear to be very sensitive to its value. The RMF-PC forces NL-Z2 [37] and 
NL3 [38], which deliver good (and also comparable) results for finite nuclei 
have quite different values, namely K = 172 MeV for NL-Z2 and K = 260 
MeV for NL3. Calculations of excited states, however, favor the value of NL3. 
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5 Conclusions 



We have presented an extended NJL model for both nuclear matter and finite 
nuclei. We tested the compatibility of the ansatz and the NJL mechanism of 
mass generation due to spontaneous breaking of chiral symmetry with nu- 
clear bulk properties and ground-state observables of finite nuclei. We found 
compatibility in both cases. 

In nuclear matter, we investigated the effect of a small (160 MeV) explicit 
chiral symmetry breaking mass term as well as the cutoff dependence. In least- 
squares adjustments, a cutoff of A c ~ 2.0 — 2.75 fm _1 is favoured compared 
to higher values. It remains to be understood why the nucleon NJL model 
favours such small values. Its similarity with the a meson mass value required 
in RMF-FR models deserves further attention. An explicit mass term of 160 
MeV leads to better results compared to the chiral limit of the model. This 
corresponds to the physically realistic case. A too small cutoff and a large 
explicit mass term do not appear to be consistent with dynamically generated 
mass. 

The results for finite nuclei show that a simultaneous description of the nucleon 
mass and ground-state observables is possible: reasonable values of mo and A c 
could be found. However, the overall quality of agreement with the measured 
observables utilized in the least-squares adjustment is not as good as that 
obtained in the RMF models. 

An explicit mass term, the additional interaction terms, and the adjustment 
procedure with respect to finite nuclei enhance the performance of the model 
to a great extent. The fact that a not too large cutoff and m = 160 MeV are 
favorable can be understood from the experience with RMF-PC models: in 
the limit m = 938.9 MeV and A c = fm" 1 we obtain the RMF-PC model, 
apart from pairing, with one exception: a 6 fermion term is not present due 
to the demands of chiral invariance of the interaction terms of the NNJL 
lagrangian that we present. The NNJL model in all cases tends to a rather 
high compressibility 310—390 MeV) in nuclear matter. Though this value is 
too high, the description of finite nuclei is still fairly good. The compressibility 
problem requires further study. 

As expected, the model for finite nuclei with a small effective mass m* per- 
forms better than the model with a larger effective mass: accurate spectral 
properties demand large scalar and vector potentials and hence a small effec- 
tive mass. Since the mean-field potential is only sensitive to the sum of scalar 
and vector potentials, this situation could be changed with terms that modify 
the spin-orbit potential, e.g. tensor terms. For the sets adjusted with respect 
to finite nuclei, the violation of causality at high densities prohibits their ap- 
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plication far above saturation density, and we attribute this behavior to our 
phenomenological approach. 

In its present form the NNJL model performs reasonably, with the too large 
effective mass and compressibility perhaps being its major drawbacks. How- 
ever, there does not appear to be much room for improvement to remedy these 
drawbacks. Obvious refinements include (a) a smooth 3-D or a covariant 4-D 
cutoff for the negative energy states, (b) pairing, (c) extended adjustment pro- 
tocols, and (d) additional interaction terms (including different choices of m , 
for example the pion-nucleon sigma term) , but it is not clear that the size and 
density dependence of the potentials will significantly change. In contrast to 
qNJL models, the applicability of this model appears to be confined to nuclear 
systems in the mean-field approximation. Nucleonic degrees of freedom and 
the current predicitve power of the NNJL model may not allow an extension 
to further areas of nonperturbative QCD. 

An objective of this study was to test if the NJL mechanism of dynamic mass 
generation due to spontaneous symmetry breaking of linearly realized chiral 
symmetry translates, in some way, from quarks to nucleons. Our correspnding 
NNJL model has not failed in this endeavor, but instead has had limited 
success. This suggests that pursuit of the linkage between low-momentum 
QCD and the nuclear many-body problem should be continued. However, 
with regard to the description of nuclear ground-state properties and the use of 
QCD scales, the nonlinear realization of chiral symmetry is superior compared 
to linear models like the present NNJL approach or the existing a — uo models. 
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